Large scale fluctuations and dynamics of the Bullard - von 

Karman dynamo 

Gautier Verhille, Nicolas Plihon, Gregory Fanjat, Romain 
Volk, Mickael Bourgoin, and Jean-Frangois Pinton 
Laboratoire de Physique de I'Ecole Normale Superieure de Lyon, 
CNRS & Universite de Lyon, F-69364 Lyon, France 

Abstract 

A synthetic fluid dynamo built in the spirit of the Bullard device |l| is investigated. It is a 
two-step dynamo in which one process stems from the fluid turbulence, while the other part is 
an alpha effect achieved by a linear amplification of currents in external coils [3]. Modifications 
in the forcing are introduced in order to change the dynamics of the flow, and hence the dynamo 
behavior. Some features, such as on-off intermittency at onset of dynamo action, are very robust. 
Large scales fluctuations have a significant impact on the resulting dynamo, in particular in the 
observation of magnetic field reversals. 
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I. INTRODUCTION 



a — uj dynamo models are extensively used in the description of natural dynamos . In 
these approaches, the u- effect converts a poloidal field into a toroidal component, while 
the a-effect is invoked for the conversion back from toroidal to poloidal. In planets, the 
electrically conducting fluid is a molten metal (iron for the Earth), i.e. a fluid with a very 
low magnetic Prantl number. As a consequence, turbulence has to be fully developed for 
induction effects to be able to overcome (Joule) dissipation, and hence induce a dynamo. 
Turbulent dynamos have been extensively studied in the context of mean-field hydrody- 
namics {4, 5]. In this approach, the turbulent small scale induction processes contribute 
to the dynamics of the large scale magnetic field through the mean electromotive force 
Zi = aijBj + PijkdjBk (the reader is referred to [6| for generalised expressions). The a and j3 
tensors depend on the small scale turbulence, and expressions for them have been established 
in several contexts. Measurements, however, have revealed that their contribution is small, 
at least at low mag netic Prandtl nu mb e rS and derate mag net 1C Reyno.ds nu.be, BQ. 
An a contribution can also originate from large scale motions through Parker's stretch and 



twist mechanism 



0,Q. 



Recently, the observation of a self-sustained axisymmetric 
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14, 15]. 



netic field in the VKS experiment [ll| has been modeled as an a-u dynamo 
This experiment has also revealed a rich time-dynamics of the magnetic field [16|, 
raising issues concerning the contribution of flow fluctuations on the dynamo behavior. It 
has been shown indeed in several studies that von Karman flows display large fluctuations 
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211 ] which have a leading role on magnetic induction 



at very low frequency in time 
characteristics 22]. 

In the present study we consider experimental synthetic dynamos 2J in which an a-effect 
is modeled by means of an external wiring, while a second induction process originates from 
the flow itself and takes into account the full turbulent fluctuations. One of these synthetic 
dynamos is an a — u dynamo similar to the one described in Q]. The u effect is linked to the 



flow differential rotation and has been extensively studied in previous experimental 22|, [23] 
and numerical works {9J : if a magnetic field is applied parallel to the cylinder axis by a set 
of coils, the fluid differential rotation induces a toroidal magnetic field. Since the flow is fully 
turbulent (with integral kinetic Reynolds exceeding 10 5 ), the induced toroidal magnetic field 
is highly fluctuating: the u-e&ect incorporates the flow turbulence. An effective a-effect is 
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then achieved when the intensity of the induced toroidal field is used to control a linear 
amplifier driving the current in the coils. This feed-back loop generates a dynamo. 

The main goal of the present study is to characterize how the resulting dynamo bifurcation 
and magnetic field dynamics are affected when large scale properties of the von Karman flow 
are modified by the insertion of various appendices inside the vessel 24J, or when rotating 
the impellers at unequal rates - as a way to impose global rotation onto the flow 191 ] . 

The experimental set-up is described in details in section [Til The dynamics of the syn- 
thetic Bullard-von Karman dynamos are described and analyzed in section III. Conclusions 
are drawn in section IV. 



II. EXPERIMENTAL SETUP 
A. von Karman gallium flow 

The synthetic Bullard-von Karman dynamo is built upon a von Karman flow. This flow 
is produced by the rotation of two impellers inside a stainless steel cylindrical vessel filled 
with liquid gallium. The cylinder radius R is 97 mm and its length is 323 mm. The impellers 
have a diameter equal to 165 mm and are fitted to a set of 8 blades with height 10 mm. 
They are separated by a distance H = 203 mm. The impellers are driven by two AC-motors 
which provide a constant rotation rate in the interval (i 7 !,^) G [0.5,25] Hz. In most of 
the cases, the flow is driven by symmetric forcing at i*\ = F 2 = F. The system is cooled 
by a water circulation located behind the driving impellers; the experiments are run at a 
temperature between 40°C and 48°C Liquid gallium has density p = 6.09 x 10 3 kgm -3 , 
electrical conductivity a = 3.68 x 10 6 ohm _1 m _1 , hence a magnetic diffusivity A = l//ioO" = 
0.29 m 2 s _1 . Its kinematic viscosity is v = 3.1 x 10~ 7 m 2 s _1 . The integral kinematic and 
magnetic Reynolds numbers are defined as Re = 2nR 2 F/u and R m = 2irR 2 F/\. Values of 
R m up to 5 are achieved, with corresponding Re in excess of 10 6 . 

Magnetic induction measurements are performed using Hall sensor probes inserted into 
the flow in the mid plane. Data are recorded using a National Instrument PXI-4472 digitizer 
at a rate of 1000 Hz with a 23 bits resolution. 
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FIG. 1: Experimental setup : a von Karman flow of liquid gallium is generated in a cylindrical 
vessel between two counter rotating impellers driven by two AC-motors. (a) Synthetic a — uj 
dynamo (axial dipole), the turbulent induced azimutal field is measured at 0.56 R in the mid-plane 
between impellers, (b) Synthetic a — BC dynamo (equatorial dipole), the turbulent induced axial 
field is measured at 0.87 R in the mid-plane between impellers, (c) Cylinder (C) configuration of 
the vessel and coordinate axis, (d) Ring {R) configuration and (e) Baffles {B) configuration. 

B. Dynamo loops 



Axial dipole 

A first dynamo configuration which generates an axial dipole is shown in figured] (a). As 
the flow is forced by counter-rotating impellers, the differential rotation in the fluid motion 
induces an azimutal magnetic field Bg when an applied axial magnetic field is imposed 
by axial coils; this is the cj-effect {4]. This induced azimuthal field is measured from a local 
probe in the mid plane at r = 0.56-R and its value is used to drive the current source which 
feeds currents into the axial coils. This external amplification constitutes the a part of the 
dynamo cycle, since an axial field stems from an azimuthal one. The feed-back amplification 
is linear (icons oc B\ up to a saturation value /foils' which is the maximum current that can 
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be drawn from the current source). The main component of the resulting dynamo is an axial 
dipole. 

For all available values of the rotation rate of the driving impellers, the magnetic 
Reynolds number is low enough to be in a quasi-static approximation. In this regime, 
Bq ~ (A^ 1 / \)B A d z vg, where vg is the azimutal velocity [9|. The u part of the dynamo 
cycle is thus very much dependent on the statistical properties of the velocity gradients at 
all scales. 



Equatorial dipole 

Another possibility is to arrange the external coils so that they generate a field perpendicular 
to the rotation axis, say B x (an equatorial dipole). In this case the induced magnetic field 
components which drive the current source is the axial induction B^J 2 measured in the 
mid plane at 90° from the applied field. Numerical studies 25| have shown that, in the 
quasi-static approximation, B?J 2 is mainly induced from B x and the radial vorticity u r of 
the flow via a boundary-condition effect due to the jump in electrical conductivity at the 

wall: Bl 12 (Vx)- 1 (A- 1 /A)([A] J B .V)c; r , where [A] is the jump in magnetic diffusivity 

at the vessel boundary. In the following we will refer to this mechanism as .BC-effect, for 
Boundary- Conditions effect, and the equatorial dipole dynamo is of a — BC type. 

In von Karman flows, the .BC-effect is particularly intense because large scale radial 
vortices develop in the mid plane as a result of a Kelvin-Helmholtz like instability in the 
shear layer 20, [2JJ] . These radial vortices are moreover known to have a highly fluctuating 
dynamics which in turn generates strong fluctuations of the BC induced field. 



Implementation 

For each configuration, the gain G, of the linear amplifier which imposes the current output 
in the coils from the induced magnetic field value is kept constant. The control parameter 
being the rotation frequency of the impellers Fi, we adjust the value of G so that the critical 
rotation frequency for dynamo onset is F c = 8.3 Hz. This gain is based on time-averaged 
values of the magnetic induction when a fixed current is imposed in the external coils - an 
'open-loop', or induction configuration. One first measures the efficiency A e g of induction at 
the probe location by (-B 7 )t = A c ffF(B A ) t , where (.) t is the time-averaged operator. Then 
one sets the gain G of the closed-loop conversion (dynamo configuration), B A = GB 1 such 
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that the estimated dynamo threshold F c = l/GA e s is constant and equal to 8.3 Hz. In other 
words, the gain is adjusted such that, for a given flow forcing parameter F, the time-averaged 
value of the amplified induced field G(-B 7 ) t is constant whatever the flow configuration and 
dynamo type. 

Saturation is set by the maximum current available from the source; the saturated mag- 
netic field value is of the order of a few tens of Gauss. As a result, the Lorentz force is too 
weak to react back onto the flow field (the interaction parameter is less than 10 -2 ). The 
synthetic Bullard-von Karman is thus similar to a kinematic dynamo. 



C. Flow configurations 

The scope of the present study is to investigate the influence of the flow dynamics on the 
behavior of the self-sustained magnetic field. Modifications are achieved by changing the 
flow, either by inserting appendices to the inner wall of the cylinder or by driving the flow 
assymetrically (i.e. rotating the impellers at different speeds). 

Three geometric configurations have been tested. They are shown in figure [T^c,d,e). 
In the base configuration (hereafter called C for Cylinder) there is no attachement. In 
configuration B, four longitudinal baffles are attached; they are rectangular pieces 150 mm 
long, 10 mm tall and 10 mm thick, mounted parallel to the cylinder axis. In configuration R, 
a ring is attached in the mid-plane; it is 15 mm tall and 4 mm thick. The respective influence 



of these attachments have been studied in water experiments [24| where it was found that 
they modify significantly the dynamics of the differential shear layer in the mid-plane. The 
presence of the ring or the baffles tends to stabilize the central shear layer; they modify also 
the number and time evolution of the radial vortices in the shear layer. 

Similarly, when the flow is driven by counter-rotating the impellers at different speeds, the 



flow bifurcates between 2-cells and 1-cell configurations |l9l.l2l|. hence leading to substantial 
modifications of the flow turbulence and large scale topology. The unbalance in the driving 
of the flow is measured by the ratio 9 = (Fx — F 2 )/(F 1 + F 2 ), 6 = corresponding to 
exact counter-rotation and = 1 to the fluid being set into motion by the rotation of one 
impeller only. The transition occurs at a critical value of above which the shear layer 
abruptly moves from the mid-plane to the vicinity of the slower impeller (the flow transits 
then from a two counter-rotating cells geometry to a single cell geometry). It is associated 
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with an abrupt change in the mechanical power required to keep the impellers rotating at 
fixed rates. 

Altogether, by inserting appendices in the flow or driving it asymmetrically, one has 
the possibility to vary the mean profiles and fluctuation characteristics of the von Karman 
flow. The impact of these changes on the synthetic a — u and a-BC dynamos is studied in 
detail below. Note that though both u and BC effects are primarily generated by the same 
differential rotation, they are sensitive to different topological properties, with also different 
dynamical behavior of the mid-plane shear layer: the 5C-effect is sensitive to the radial 
vorticity (u r ) in the vicinity of the vessel wall while the u;-effect is sensitive to the axial 
change of the toroidal velocity in the bulk of the flow. The natural fluctuations of these flow 
gradients are known to be of different kinds : radial vorticity fluctuations are closely linked 
to the dynamic of strong intermittent Kelvin-Helmhotz vortices in the shear layer, while the 
axial gradient of toroidal velocity fluctuates according to local and global displacements of 
the shear layer. Exploring both the axial and equatorial dipoles configuration give us then 
a way to probe selectively the impact of these velocity gradients on the dynamo. 

III. DYNAMO AND FLOW DYNAMICS 
A. Axial dipole, a — u dynamo 

In this section, the flow is driven symmetrically: impellers counter rotate at the same rate 
F = F 1 = F 2 . The external coils are set to generate an axial dipole field. The turbulent uo- 
effect (the source of the fluctuations for the dynamo) is first analyzed as an induction process 
on its own. Dynamical features of the a — u dynamos are then described and analyzed. 

1. Turbulent u induction 

We first consider induction measurements obtained when a constant axial magnetic field 
is applied. We show the temporal dynamics of the induced field in figure EJ^b), and display 
the probability density functions (PDFs) and time spectra of Bg/ (Bg) in figure [2]^c-d). Note 
that in order to compare the fluctuations for the different configurations, we have normalized 
the induced field by its mean amplitude in figure [2](c,d). Indeed, in the dynamo loop, the 



7 



a) 




y 



0.2 
!o.1 



0.2 

8 

Lai 

I 

0, 



5 i r> 15 20 25 30 
-i 1 1 1 1 



10 15 20 25 30 



R' 



b) 



10 15 20 25 30 

Time (s) 



( 



10 
10° 



■ Cylinder 
« Baffles 

■ Ring 


- i ' ! 1 ! M 1 ' ! 1 


■ 1 1 ! 1 1 ! F 




Hl[jl!Hf 






■|lt|l|H/l!IW.*!!l!llt!>ll 

: : : ! : > ■ w '■ : \ : °i \ : : -: ! : « 1 t! 


j%h^| j .1 i | t e r 


)f 

■ii i m l1 1 1 1 1 1 h ! n ] 1 1 1° n ft 1 1 g 1 1 1 : r 




i ** 


■| 1 1 il i 1 1 1 1 h i n n j 1 1 1 f i f i 


ft i \ii n i r 







B « / < B B> 





Cylinder 




Baffles 




— -Ring 




FIG. 2: Axial applied field, azimutal induction features at r = 0.56 R. (a) Experimental setup, 
(b) Temporal dynamics for the three flow configurations, (c) Probability density functions and (d) 
Power spectrum density of the azimutal induced field (log-lin representation in the inset). Note 
that, for the last two graphs, the induced fields are normalized by their time-averaged values. 

value of the gain G is adjusted so that the time-averaged amplified induced field for a fixed 
rotation rate is the same whatever the flow configuration. 

• C configuration. The shear layer is fully developed in this less constrained state. As 
a result, its motion explores a large portion of the flow volume. These motions result 
in large variations of the amplitude of the induced field, with a quasi- Gaussian statis- 
tics. The power spectrum of fluctuations in time reveals an important low frequencies 



dynamics. These features have been discussed in detail in previous studies 



22] 



R configuration. When a ring is installed in the mid-plane, the shear layer is 
pinned 24j]. The axial gradient of the time-averaged azimutal velocity is peaked at the 
ring position ^24\ and higher than in the C configuration. The time-averaged azimutal 
induced field is nevertheless of the same order of magnitude in both configurations 
(cf. figure E]^b). This is because the magnetic field 'filters' the velocity gradients at a 
length scale of the order of the magnetic diffusion length r\B ~ a/A/F 5 cm at R m = 2. 
The magnetic field is more sensitive to the difference of velocity on either side of the 
shear layer rather than in the actual slope of the velocity profile. 
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Another effect is that fluctuations in the azimuthal induction are reduced due to the 
shear pinning. The time spectra for the two configurations are similar in the inertial 
range, but the presence of the ring reduces low frequency fluctuations. The slight peak 
at very low frequency is attributed to a 'depinning' occuring at irregular time intervals, 
also observed as a bimodality in the probability density function in figure [21(c) . 

B configuration. When baffles are inserted in the vessel, vortices of the shear layer 
are confined to inter-baffles position and fluctuations of the shear layer are also 



damped 



241 ] . Like in the R configuration, the time-averaged induced field is weakly 



modified, but fluctuations are highly reduced as compared to the C configuration. 



For these three configurations, the mean induced field and normalized fluctuations for 
R m = 2 are the following: (Bg/Bf) = 0.08; 0.12; 0.07 and (B^) imB / (Bg) = 0.32; 0.11; 0.1 
for the C, B and R configurations respectively. Comparison of the spectra [2]^d) show that 
modifications of the flow strongly modify the low-frequency, large scale dynamics of the 
turbulent induction, while keeping small-scale (inertial range) dynamics almost unchanged. 

2. Dynamo: axial dipole 

We then analyse the behavior of the synthetic a — uj dynamo [i.e. when the measured 
azimuthal field actually drives the current in the induction coils). Time signals of the 
axial magnetic field (normalized to the maximum field) are shown in figure [3] for the three 
configurations and an increasing forcing parameter F (from top to bottom). For the three 
configurations and whatever the value of the forcing parameter, the axial dipole dynamo 
is always observed homopolar: when growing in either polarity, the field remained of same 
polarity. However both polarities were observed for different realizations {i.e. when stopping 
the motors and increasing the rotation rate above the critical value), as is expected from 
the B — > — B symmetry of the MHD equations. 

Let us first analyze the C configuration, for which last section showed that fluctuations 
of the u-e&ect are the highest. At very low rotation rates F <C F c , the magnetic field is null. 
When the rotation rate increases and approaches the threshold (F < F c ) , the magnetic field 
remains null most of the time, but exhibits intermittent bursts. The density and duration 
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FIG. 3: Axial dipole dynamo temporal features: temporal evolution of the dipole amplitude as a 
function of driving frequency for the three configurations. 

of the dynamo state increase rapidly with the impellers rotation rate: as F > 1.3F C , the 
magnetic field in mostly non-zero, with intermittent extinctions for short time intervals. The 
system thus bifurcates from a non dynamo to a dynamo state via an on-off intermittency 



2fl 



27|. This kind of 



scenario, also previously observed in numerical models of the dynamo 
bifurcation is characterized by specific features of the signal: its probability density function 
is peaked at zero near onset and then decays algebraically; the PDFs of 'off' times scales as 



-3/2 



off 



28]. 



Evolution of the measured PDFs of magnetic field as a function of the forcing parameter 
is shown in figure S](a). The PDF is peaked at zero for rotation frequencies up to 1.2F C , and 
extends to large values, even at low rotation frequencies. A log-log representation is shown 
in inset of figure H](a) and shows a typical negative power- law for the PDF (with an exponent 
proportional to the control parameter). This feature was observed in numerical simulation 



of on-off models 



29 



30| and simulations of MHD equations 26] . On-off intermittent signals 
also display an universal behavior of inter-burst time intervals, or time interval between 
dynamo activity, named T g - for 'off' phases or laminar phases. This time interval is defined 
in the inset of figure H^b): one sets a threshold and computes the time T Q g for which the 
dynamo field is below this threshold. The probability density functions of T Q g are displayed 
in figure H](b) for four values of the forcing parameter. These probability density functions 
have a power law dependence with a cut-off at large values (the distribution of the 'off' 
phases has been checked to be independent of the value of the threshold). A T off 3 ^ 2 scaling 
is shown for comparison with the experimental data, showing that the data statistics is in 
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agreement with the -3/2 power law scaling behavior predicted for on-off intermittency 281] . 



observed in numerical simulations 



26 



29 



301 ] or in experimental investigations in electronic 



circuits 3l|, gas-discharge plasmas 32j or convection in liquid crystals 33j. Corrections 



to the power-law scaling are not obvious, so that more complex scenarii (such as ' in-out' 



intermittency [34|]) may not apply here. 

The variation of the cut-off time in the PDF of T g- with the control parameter is steep. 
We observe a power law with a —5 exponent in a small range F G [0.97-F c , 1.69-FJ. This 
is in contrast with numerical models with multiplicative noise for which a decay as F~ 2 is 
quoted [29] - 
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FIG. 4: Axial dipole dynamo for the C configuration: (a) evolution of the probability density 
function of the axial magnetic field as a function of the driving parameter - inset: log-log plot 
for F = 0.84.F C - and (b) - inset: definition of inter-burst time interval - main: evolution of the 
probability density function of the inter-burst intervals as a function o the driving parameter. 

Modifications of the u;-effect with the other two configurations drastically change the 
dynamo features. As can be seen in figure [3] for the B and R configurations, no dynamo 
bursts are observed for impellers rotation rates well below threshold and 'off' phases are 
no longer present well above threshold. The associated bifurcation curves are displayed 
in figure 0(a). The bifurcation is steeper for the B and R configurations as compared to 
the C configuration, i.e. the width of the intermittent domain decreases when inserting 
appendices in the vessel. The on-off intermittent regime is nevertheless always observed in 
a very narrow range near onset. Figure [5^b) shows the probability density functions of the 
'off' phases for the three configurations near the critical rotation rate (F = 1.08F C ). The 
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-3/2 power law scaling is independent of the configuration, while the exponential cut off 
increases when the low frequency component of fluctuations decreases. The mean length of 
'on' phases also significantly increases when low frequency of the noise decreases, as can be 
observed in figure [3J 
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FIG. 5: Features of the axial dipole dynamo: (a) evolution of the time-averaged axial magnetic 
field as a function of the rotation frequency for the three configurations; (b) probability density 
function of the inter-burst time interval for the three configurations at F = 1.08-F C . 



As a partial conclusion, we observe that several features (such as the on-off scenario and 
waiting times statistics) of the axial dynamo at onset are robust with respect to large scale 
modifications of the von Karman flow. Other features (such as the width of the transition 
and the frequency of incursions to or from dynamo states) are on the contrary strongly 
dependent on the large scale dynamics of the flow. 



B. Equatorial dipole, a — BC dynamo 

We now consider the a-BC dynamo configuration, generating an equatorial dipole field, 
perpendicular to the axis of rotation of the impellers. As for the axial case, we study the 
dynamo bifurcation in the C, B and R configurations, the flow being driven symmetrically 
with the impellers counter rotating at the same rate F = F\ = F?,. As for the previous 
section, the gain of the feedback loop is adjusted for each configuration in order to keep 
constant the dynamo threshold estimated from time-averaged induction measurements (F c = 
8.3 Hz). 
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1. Turbulent BC induction 



Induction measurements are performed with a constant field, applied perpendicularly to 
the cylinder axis. The induced field in the axial direction is measured in the mid plane at 
90° from the applied field, at distance 0.87i? from the axis. 

• C configuration. The induced field (figure HD(b)) shows a highly fluctuating dynamics. 
As seen in figure E](c), the induced field distribution is bimodal, with a peak near 
zero induction and one at about twice the time-averaged value. The induction follows 
the dynamics of the radial vortices formed in the shear layer. A significant induced 
field is measured only when radial vortices sweep the probe. In between these events, 
the induced field is weak with large fluctuations - the sign of the induced field is not 
prescribed. 

• B configuration. As compared to the above situation, the induced field is always non 
zero, a feature consistent with the observation 24| that the radial vortices are now 



confined between the baffles. A concurrent feature is that the fluctuations of induction 
are much reduced (a factor of 5 lower than the C configuration). The spectra displayed 
in figure [6] show that this reduction occurs at all scales. 

R configuration. In this case, the shear layer is pinned. The induction signal is again 



stationary. This is consistent again with observations [2J, [35J that the vortices are 
attached to the central ring. The fluctuations are reduced again by a factor of 5, as 
also seen in figure EJ^c), this reduction lies essentially in the low frequency motions 
(figure EJ^d)). 

For the three configurations, mean induced field and normalized fluctuations at R m = 2 
are the following: (B*/B£) = 0.029; 0.029; 0.058 and [B\ ) rma /(^ ) = 0.90; 0.17; 0.03 for the 
C, B and R configurations respectively. 



2. Dynamo: equatorial dipole 

Time recordings of the equatorial dipole dynamo field are shown in figure E(a). In each 
configuration, time signals are shown for three values of the rotation rates of the impeller, 
from below the critical frequency to above. 
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FIG. 6: Transverse induction, (a) experimental set-up, (b) temporal evolution of axial induction 
Bl(t) for a constant transverse applied field (note the difference in y-axis scales) (c) probability 
density function and (d) time spectra of normalized induced fields - log-lin representation in the 
inset. 

When compared to the previous axial dynamo cases, several features emerge. First in 
the absence of appendices in the vessel (configuration C), a dynamo field grows with either 
polarity and reverses spontaneously. Below threshold, bursts with both polarities ±B cq are 
observed. For higher values of the rotation rate, the fraction of time spent in a dynamo state 
increases. This is in sharp contrast with the B and R configurations for which homopolar 
dynamos are always generated after a rather abrupt bifurcation. In these configurations, 
bursts of dynamo action are observed near threshold but a steady state with a definite 
polarity is rapidly reached above threshold. 

Let us describe and analyze in details the C configuration for which the magnetic field 
exhibits reversals. The bifurcation curve, displayed in figure [7J(b)-black shows a smooth 
transition from the null field regime to the dynamo regime. This transition occurs through 
bursts of both polarities below and around threshold. The time spent in the dynamo regime 
increases with the forcing parameter. As compared to all other Bullard-von Karman dy- 
namos obtained so far, a striking feature is that the reversing dynamo displays significant 
'off' phases even well above threshold. The probability density functions of the 'off' phases 
for three values of the forcing parameter are shown in figure [7(c). The observed statistic is 
consistent with a Poisson distribution having a characteristic time scaling as {F/F^" 1 . On 
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the other hand, the distributions of the duration of the saturated ('on'; plus or minus) states 
(not shown) are consistent with a sum of two Poisson distributions with characteristic times 
in a ratio of 7 and scaling as Fj F c . Note that in this configuration, statistics are converged 
only when working with very large time recordings. Regimes at large F/F c have thus been 
obtained by increasing the gain G in order to keep F in the optimal range for operation 
with the liquid gallium kept at constant temperature (around 40° C). 

For the B and R configurations, the dynamics of the homopolar dynamo is very similar to 
the previous axial case. The on-off intermittent behavior is controlled by the low frequency 
content of the induction process: it is here restricted to a very narrow range of F around 
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FIG. 7: Equatorial dipole dynamo behavior: (a)time signals for increasing values of the rotation 
rate of the impellers, from below to above threshold, (b) bifurcation curves and (c) probability 
density function of the 'off' phases for the C configuration, i.e. the reversing dynamo. 



As a partial conclusion here, the equatorial dynamo also develops via an on-off intermit- 
tent regime, but some features have been significantly changed, such as the ability of the 
dynamo to undergo reversals of polarity. These reversals are only observed when strong 
fluctuations are present. This is a common feature of reversals dynamics when caused by 



stochastic transitions between the symmetric B and — B states 



36] 



15 



C. Global rotation 



In this section, we adress the case where the flow is driven with impellers counter-rotating 
at different rotations rates in the presence of a ring. This is an asymmetric R configuration, 
as for the dynamical regimes of the VKS dynamo reported in Q, Q, Q. 

In order to 

emphasize the influence of the dynamics of the mid-plane shear layer, and its links with 
reversals, we focus on the equatorial dynamo loop. 



1. BC induction with asymmetric forcing 

As = (Fi — F 2 )/(Fi + F 2 ) is varied, one observes a clear transition in the torques 
Tj driving the two impellers. Figure [8](b) shows the dimensionless differences between the 
torques of the two impellers (these torques are deduced from the currents measured by the 
electrical drives that feed the motors). The bifurcation between a 2-cell flow and a 1-cell 
flow occurs at a criticl value C ~ 0.16, in agreement with measurements made in water 
flows 
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371 ] . For |0| < C the time-averaged flow consists of two main cells on either 
side of the mid plane, in which the toroidal and poloidal flows have opposite directions. 
Note that the regimes previously described correspond to a symmetric forcing, i.e. = 0. 
For |0| > C the flow volume is dominated by one cell driven by the fast impeller, as is 



schematically drawn in figure[H](b). It was shown in previous studies in water flows 19i.l24l.l37| 
that fluctuations of the flow diverge at the transition between the two flow regimes. 

Figure [8](a) shows that fluctuations of magnetic induction strongly depend on the asym- 
metry of the forcing. Energy distribution across scales of the turbulent induced field also 
strongly depend on the asymmetry of the forcing as can be seen in figure [8]^c). Near ~ 0.16 
(shear layer instability) the low frequency part of the spectrum is enhanced. At higher 
values the spectrum is enhanced around and above the forcing frequency. One thus estimates 
that the flow has slow large scales fluctuations only near 0°, while it has larger small scale 
fluctuations for 3> 0°. 

Evolutions mean induced field and normalized fluctuations at R m = 2 as a function 
of global rotation are the following: (B T g /B^) = 0.058; 0.040; 0.013 and [B\ ) nns /(B J / ) = 
0.03; 0.16; 0.22 for the = 0; 0.16; 0.6, respectively. 
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FIG. 8: Influence of global rotation on BC induction features: (a) time signals for three values 
of the global rotation, (b) dimensionless torque difference as the function of the flow asymmetry 
- insets show sketchs of the mean von Karman flows, (c) time evolution of the induced field, and 
(d), corresponding power spectra. 

2. Asymmetric equatorial dynamo 

Time evolution of the equatorial dynamo field B x is shown in figure [9] for three values of 
G (refer to figure [7] for the = case) and increasing values of the rotation rate F. 

O X TT0 V C 

= O.16j \ » = 0.2 | \ % © = 0.6 \ » 




Time (s) Time [s| Time (5} F/Fc 

FIG. 9: Transverse dynamo with global rotation: time recordings of the transverse magnetic field, 
for increasing values of G (left to right) and for increasing value of the mean rotation rate F. 
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The = case, described in last section, showed an homopolar dynamo with an on- 
off intermittent regime at onset and a steep transition from the null-field regime to the 
dynamo regime. When increasing 0, a first observation is that both polarities occur, with 
chaotic reversals for ~ C . This feature disappear for higher values of 0, i.e. when 
the flow has bifurcated toward the one-cell configuration, and one recovers an homopolar 
dynamo. The analysis of the ~ C regimes is of great interest since, for these regimes, 



the VKS dynamo displayed most o 
asymmetric bursts, oscillations...) 



its dynamical regimes (chaotic reversals, symmetric and 
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17| . At the critical value ~ C , the equatorial 



dynamo is bipolar, bifurcates through an on-off scenario and reverses chaotically. Close to 
the threshold value (the = 0.2 regime is displayed in figure ED, the dynamo also bifurcates 
via an on-off scenario and exhibits bursts of both polarities for a rotation rate close to F c . 
At higher rotation rates, reversals occurences are less frequent and the dynamo eventually 
reaches an homopolar regime for F/F c > 1.5. 

When the fow has bifurcated to the one-cell configuration, the equatorial dynamo is 
homopolar. On/off intermitency is observed in a very narrow range of driving frequencies. 



IV. CONCLUDING REMARKS 



In the present study, influence of the flow fluctuations on synthetic Bullard-von Karman 
dynamos has been investigated. Two types of dynamo loops were studied: ana-u and 
a a — BC dynamo, where the u and £>C-effects incorporate turbulent fluctuations. Flow 
fluctuations modifications have been achieved by inserting appendices in the vessel, or by 
driving the flow asymmetrically (for a — BC dynamo). Several robust features have been 
observed: 

• The bifurcation occurs via an on-off intermittent regime at onset of dynamo action. 

• The on-off intermittent regime is controlled by the low frequency part of the fluctuating 
induction process considered in the dynamo loop: the higher the low frequencies of 
the fluctuations, the wider is the occurence of on-off intermittency among the control 
parameter. 

• For all studied configurations, the system spends half ot its time in the dynamo regime 
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when driven at the critical forcing parameter defined on time-averaged induction pro- 
cesses. 

However, some features strongly depend on the exact configuration. Reversing dynamos 
have been observed only for strong fluctuations of the turbulent induction process. The low 
frequency part of the spectrum seems to play also a dominant role on the ability of the 
dynamo to reverse. 

Hence, for this synthetic dynamo, salient features such as bursts of magnetic field activity 
or reversals are controlled by the underlying hydrodynamics. A detailed comparison with 
the dynamics of stochastic differential equations in the presence of noise is underway and 
will be reported elsewhere. 
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